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ABSTRACT 

Upcoming and ongoing large area weak lensing surveys will also discover large sam¬ 
ples of galaxy clusters. Accurate and precise masses of galaxy clusters are of major 
importance for cosmology, for example, in establishing well calibrated observational 
halo mass functions for comparison with cosmological predictions. We investigate the 
level of statistical uncertainties and sources of systematic errors expected for weak 
lensing mass estimates. Future surveys that will cover large areas on the sky, such as 
Euclid or LSST and to lesser extent DES, will provide the largest weak lensing cluster 
samples with the lowest level of statistical noise regarding ensembles of galaxy clusters. 
However, the expected low level of statistical uncertainties requires us to scrutinize 
various sources of systematic errors. In particular, we investigate the bias due to clus¬ 
ter member galaxies which are erroneously treated as background source galaxies due 
to wrongly assigned photometric redshifts. We find that this effect is significant when 
referring to stacks of galaxy clusters. Finally, we study the bias due to miscentring, 
i.e., the displacement between any observationally defined cluster centre and the true 
minimum of its gravitational potential. The impact of this bias might be significant 
with respect to the statistical uncertainties. However, complementary future missions 
such as eROSITA will allow us to define stringent priors on miscentring parameters 
which will mitigate this bias significantly. 
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1 INTRODUCTION 

Galaxy clusters play an important role in testing cosmologi¬ 
cal models, for example, by confronting the observed number 
of galaxy clusters with predictions of the halo mass function 
(e.g. Tinker et al. 2008, 2010 and references therein). This 
test is especially sensitive to the values of the matter con¬ 
tent of the Universe, 12^, and the normalization of the pri¬ 
mordial power spectrum of matter density fluctuations, ag. 
Values for both parameters obtained from recent Sunyaev- 
Zel’dovich (SZ) cluster counts by Planck Collaboration et al. 
(2014a, 2015b) are in tension with other independent mea¬ 
surements (e.g. Planck Collaboration et al. 2014b, 2015a). 
To relate the observed cluster counts to predictions for the 
mass function a conditional scaling relation was used. The 
analysis in Planck Collaboration et al. (2014a) was based on 
X-ray mass proxies. The uncertainty in the absolute clus¬ 
ter mass scale remains the largest source of uncertainty in 
the Planck cluster count analyses and is quantified by the 
mass bias. New priors on this mass bias were incorporated 
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in Planck Collaboration et al. (2015b) based on small over¬ 
lapping cluster samples with masses measured by employing 
gravitational lensing (von der Linden et al. 2014; Hoekstra 
et al. 2015), the deflection of light due to mass as a conse¬ 
quence of Einstein’s equivalence principle. Although these 
improved priors do not fully lift the observed tension yet, 
scaling relations entirely based on and not only gauged by 
gravitational lensing measurements are advantageous. This 
is due to lensing masses being unaffected by the dynamical 
state of matter or its physical properties (e.g. being dark 
or baryonic) in general. One disadvantage though is that 
lensing only yields estimates for a two-dimensional surface 
mass density, but with simulations it is possible to propa¬ 
gate these reliably into three-dimensional mass estimates in 
order to compare them with results from other probes (e.g. 
Becker & Kravtsov 2011; Meneghetti et al. 2014). 

In the strong limit gravitational lensing is characterised 
by the occurrence of multiply lensed images or arcs of back¬ 
ground sources behind the cluster. Employing these lensed 
images allows a very detailed determination of the mass and 
mass profile of the cluster core (e.g. Medezinski et al. 2013, 
Kohlinger & Schmidt 2014; Bartelmann et al. 2013 and ref- 
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erences therein). In the weak limit small differential deflec¬ 
tions of background galaxies are used in a statistical sense 
to infer the mass of the cluster (e.g. Hoekstra et al. 2013 
and references therein). 

The primary source of statistical uncertainty in weak 
lensing based cluster mass estimates is shape noise because 
galaxies are not intrinsically round. For weak lensing it is 
necessary to measure shapes accurately from observed im¬ 
ages in order to derive ellipticity components which serve 
as shear estimators (cf. Section 2.3). The accumulated er¬ 
rors arising from measuring shapes eventually propagate 
into weak lensing analyses as shape noise. 

Another source of statistical uncertainty is arising from 
the fact that mass as measured from gravitational lensing 
is always weighted by the lensing kernel along the line of 
sight and projected into the plane on the sky of the deflect¬ 
ing mass. The effect of this projected fore- and background 
mass - or cosmic noise - on the accuracy of weak lensing 
masses has already been studied extensively in the past (e.g. 
Hoekstra 2001, 2003; Dodelson 2004; Hoekstra et al. 2011). 
However and in particular for the cosmological test described 
above, one is interested in a stack of clusters within a given 
mass (and redshift) range, thus, statistical uncertainties on 
properties of the stack will scale inversely with the square 
root of the total number of clusters within these bins. 

In the next decade data from a multitude of ground- 
based weak lensing surveys, for example, the Kilo-Degree 
Survey (KiDS^, de Jong et al. 2012), the Dark Energy Sur¬ 
vey (DES^, Flaugher 2005), and the Subaru Hyper Suprime- 
Cam lensing survey (HSC^) will become available and even¬ 
tually culminate in the surveys carried out by the Large 
Synoptic Survey Telescope (LSST^, Ivezic et al. 2008) and 
the spaceborne EucliF mission (Laureijs et al. 2011). Al¬ 
though the major focus of these surveys will be cosmic shear 
- the much weaker weak lensing due to cosmological large- 
scale structure they will also produce large cluster sur¬ 
veys as ’by-products’ and allow mass estimates employing 
shear measurements (cf. Sartoris et al. 2015). Since all of 
these weak lensing surveys will provide superior statistics in 
terms of the expected number of clusters to be found due 
to an orders of magnitude increase in survey area, it is im¬ 
portant to scrutinize sources of statistical and systematic 
uncertainties in the determination of cluster masses using 
weak gravitational lensing. 

In this paper, we will be especially focussing on the 
Euclid survey because it will yield one of the largest clus¬ 
ter samples with a very low level of statistical uncertainties 
based on weak lensing. This allows us to scrutinize the im¬ 
pact of major sources of uncertainties and biases, eventu¬ 
ally answering whether these large cluster surveys will be 
more affected by statistical uncertainties or systematic er¬ 
rors. Note, however, that the survey design for LSST is very 
similar to Euclid so that our results will also be applicable 
to this survey and to lesser extent to smaller surveys such 
as, for example, DES. 

Due to the large increase in survey area, we expect large 
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samples of clusters to be detected in ongoing and upcoming 
weak lensing surveys and hence the statistical uncertainties 
are expected to decrease to levels on which systematic errors 
will no longer be negligible but instead might even dominate 
over statistical uncertainties. As we have mentioned already, 
gravitational lensing always yields a projected mass due to 
the line of sight integration over the lensing kernel. Prefer¬ 
ably, the redshifts of (all) background source galaxies should 
be used in the integration over the lensing kernel (cf. Sec¬ 
tion 2.3) which thus requires to take redshifts of even larger 
samples of galaxies. In the ongoing and upcoming surveys 
this will be achieved by employing photometric redshift es¬ 
timates based on multiband observations since spectroscopy 
is not feasible anymore given the typical survey areas (at 
least several thousand square degrees). 

However, photometric redshift estimates are less precise 
and accurate and photometric misidentifications, for exam¬ 
ple, in the cluster member galaxy assignment will propagate 
into a biased mass estimate via the lensing kernel. Similarly, 
the effect of miscentring - the displacement between any 
observationally defined cluster centre and the minimum of 
the gravitational potential of the cluster - can also be prop¬ 
agated into a bias of weak lensing mass estimates. On the 
scale of a single cluster these biases are negligible but this 
might not be anymore the case once we turn to studying the 
masses of large ensembles of clusters. 

In this paper, we will determine the level of expected 
statistical uncertainties first and then continue to study pos¬ 
sibly important sources of bias. Eventually, these biases will 
have to be assessed much more rigorously through extensive 
(numerical) simulations, the aim of this paper is, however, 
to provide a guideline for the design of these simulations by 
identifying the most significant sources of bias with respect 
to the expected level of statistical uncertainties for stacks of 
galaxy clusters. 

The structure of the paper is as follows: in Section 2 
we establish the level of expected statistical uncertainties 
on the mass estimates of stacks of galaxy clusters from the 
Euclid survey while describing the weak lensing formalism 
used in this study at the same time. In Section 3 we scru¬ 
tinize various sources of bias, most importantly the effect 
of photometric redshift outliers as well as miscentring and 
compare them to the level of statistical uncertainties derived 
in the previous section. Finally, we present our conclusions 
in Section 4. 

Throughout this paper we employ a spatially flat 
ACDM cosmology with = 0.308, Qa = 0.692, Hq — 
lOOhkms”^ Mpc“^ with h = 0.678, erg = 0.826 and Ua = 
0.961 following results from Planck Collaboration et al. 
(2014b). 


2 STATISTICAL UNCERTAINTIES 

Estimating the total number of haloes in redshift and clus¬ 
ter mass bins is the first step in answering the question of 
whether a Euclid cluster survey will be limited by statistics 
or systematics since the statistical uncertainty on a stack 
of clusters scales inversely with the square root of the to¬ 
tal number of haloes in the stack. In a real data analysis 
one has to stack the clusters according to an observational 
proxy (e.g. luminosity). This will introduce an Eddington 
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bias (Eddington 1913) in the stacked quantity as a func¬ 
tion of decreasing number density. Hence, the average halo 
masses on the high mass end will be lowered. This in turn 
reduces the signal-to-noise ratio of the weak lensing mea¬ 
surement. However, this would only weaken the constraints 
on systematic errors we derive below. Therefore, we do not 
assume a proxy for the stacking but stack the clusters di¬ 
rectly in mass in order to derive the tightest constraints on 
systematic errors. 

The Euclid survey will cover an area of 15000 deg^ on 
the sky and is expected to detect 30 galaxies per square 
arcminute (Laureijs et al. 2011) for which accurate shapes 
can be determined. For the subsequent analysis we will 
assume that all clusters within this area will be detected 
down to redshifts of 2 = 1.5 spanning masses between 
6.78 X lO^®/i“^M 0 and 2.70 x lO^^/i'^M©. This assump¬ 
tion is deliberately optimistic because it leads to the tightest 
constraints on systematic errors. Note, however, that for ex¬ 
ample Gladders et al. (2007) have found their high-redshift 
sample of clusters from the Red-Sequence Cluster Survey 
(RCS) to be complete to « 88 %. 


2.1 Halo abundance 


The halo abundance can be expressed in the functional form 
(Tinker et al. 2008, 2010) 


\ t^-Tl Pm,0 r{ \ 

The function f{v) is motivated from extended Press- 
Schechter theory and can explicitly be written as (Tinker 
et al. 2008, 2010 ) 




-2<^l 2rj —^1^2/2 


( 2 ) 


where the parameters and (j) are redshift- 

dependent and have to be calibrated against numerical sim¬ 
ulations for the corresponding overdensity A = 200 with 
respect to prn (cf. Tinker et al. 2008, 2010 for explicit values 
of these parameters). 

The halo abundance is predicted as a function describ¬ 
ing the mass fraction of matter in peaks of a given height, 
u = 5c/cy{M,z), in the linear density held smoothed at 
a scale R = (3M/(47rpm,o))^^^ (Press & Schechter 1974), 
where M and 2 refer to the cluster mass and redshift, re¬ 
spectively, and pmfi is the mean matter density of the Uni¬ 
verse today. The constant 5^ = 1.686 denotes the critical 
overdensity for collapse in linear theory and cr(M, z) is the 
root-mean-square (rms) variance of the linear density held 
smoothed on a scale R{M), which is dehned as 

P{k)\W{k- R)\'^k'^dk. (3) 

Here, P{k) is the linear matter power spectrum, which we 
calculate with the htting formulas provided by Eisenstein & 
Hu (1999) and W is the Fourier transform of the top-hat 
hlter with radius R in real space. 

Using equation 1 we can now predict the (expected) 
abundance of haloes per redshift bin i and mass bin j by 
evaluating 

i cA'Cigh j r 

Nii= / / n(M, ,2i)d2dMdUcom, (4) 



where we integrate over the expected comoving volume of 
the Euclid survey. The subscripts ’low’ and ’high’ refer to 
the lower and upper bounds of the bin, respectively. 


2.2 Mass model 


Next we have to specify a mass model for galaxy clusters 
from which we will derive shears that can then be compared 
to the measured shear around galaxy clusters. Numerical 
simulations of cosmological volumes show that the Navarro- 
Frenk-White (NFW)-proHle (Navarro, Frenk & White 1997; 
Navarro et al. 2010) is a good description of the average 
density profile of an ensemble of haloes over several orders 
of magnitude in mass when adjusting the halo concentration 
accordingly. 

In the following analysis we assume only a single halo 
component and a spherically symmetric distribution of the 
cluster mass. In general, this assumption is over-simplifying 
and especially for unrelaxed single haloes far from correct 
(e.g. Shaw et al. 2006). However, since we focus in our anal¬ 
ysis on a stacked signal from an ensemble of clusters, this 
simplification holds, because non-spherical symmetric clus¬ 
ter geometries will average out in the stacking process pro¬ 
vided that the selection of the cluster sample is unbiased. An 
unbiased cluster sample is an important assumption here in 
order to derive upper limits on systematic errors, but in a 
real data analysis the cluster selection function has to be 
fully taken into account, e.g., in a subsequent cosmological 
analysis. 

In this case the radial proHle of such an idealized halo 
can then be expressed as an NFW-profile: 


p(r) 


5 ■ pm{z) 

r/r^l -br/r ,)2 ’ 


(5) 


where pm(z) is the mean matter density of the Universe at 
the redshift 2 of the halo. The parameter S describes the 
overdensity of the halo and is related to the concentration 
parameter c through 


200 _ ^ _ 

3 ln(l -be) — c/(l + c) 


( 6 ) 


The scale radius is a characteristic radius of a cluster and 
can be related to the virial radius r 2 oo and concentration 
parameter c via = r 2 oo/c. We define the virial radius here 
as the radius of a sphere which contains a mass overdensity 
of 200pm (z). Thus, the corresponding mass M 200 within this 
sphere is given by 

M 200 = -^pm{z)r2oo ■ (7) 

Furthermore, numerical simulations hint at a (noisy) rela¬ 
tion between halo concentration and mass. By applying such 
a concentration-mass-relation, we can reduce the free param¬ 
eters of the model to only one: the mass M 200 . 

For the concentration-mass relation we use the results 
of Dutton & Maccio (2014), i.e., 

logio(c20o(.Af20o)) = n + felog2o(.^20o/(10^^ h ^AI©)) (8) 


with the redshift-dependent functions a = 0.520 + (0.905 — 
0.520) exp(-0.6172^-^^) and b = -0.101 + 0 . 0262 . The con¬ 
centration C200 and mass M 200 are defined with respect to 
the critical density of the Universe. We convert between 
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this definition and our definition of mass and concentration 
given with respect to the mean matter density employing 
the algorithm from Hu & Kravtsov (2003). Note that at 
an earlier stage of the subsequent analysis we employed the 
concentration-mass relation from Duffy et al. (2008) which 
qualitatively did not affect any of our subsequent results or 
conclusions. That is expected because the weak lensing sig¬ 
nal depends to first order on mass only. Moreover, we do 
not assume any scatter in the concentration-mass relation, 
because scatter will mainly affect the shape of the profile at 
small scales. However, our analysis always assumes that we 
measure the weak lensing signal of a stack and that these 
these small scale fluctuations from halo to halo due to scat¬ 
ter in the concentration-mass relation average out. 


2.3 Weak lensing formalism 


Analytical formulas for the calculation of the weak lensing 
convergence and shear signal from a spherically symmetric 
NFW-profile were derived in Bartelmann (1996) and are 
conveniently re-expressed in Wright & Brainerd (2000). Fol¬ 
lowing these references, we write the convergence as 


knfw(2:) 


r{x) 


(9) 


which is thus the ratio of the surface density ENFw(a;) at 
projected position x = R/vs scaled by the critical surface 
density 


E 


crit 


4-170 Di{zi) 




( 10 ) 


where = Ds{z)/Dis{z, zi) is the inverse of the lensing 

efficiency I3{z). Here Ds, Di, and Dis denote the angular di¬ 
ameter distances between observer and source, observer and 
lens, and lens and source, respectively. The constants c and 
G are the speed of light and gravitational constant, respec¬ 
tively. For explicit formulas of the surface density Enfw( 2:) 
of an NFW-profile we refer the reader to the original liter¬ 
ature (Bartelmann 1996; Wright & Brainerd 2000). 

As we have hinted at already in Section 1, weak lensing 
requires the knowledge of the redshift of the lens and every 
background source which are entering as variables in the 
corresponding angular diameter distances. However, instead 
of considering a redshift for every single background galaxy 
one assumes or measures a source redshift distribution: 


Psrc (^) 


P 

^oF(i±^) 



exp(-(2/zo)^), 


( 11 ) 


where we have adopted the functional form presented in 
Vafaei et al. (2010) and use a = 0.96, (3 = 1.70, and 
zq = 1.07 corresponding to a median redshift of ^med = 0.91 
to simulate the Euelid survey. Employing this source red¬ 
shift distribution lets us rewrite the critical surface density 
as 

= 4nGDiizi) ( 12 ) 

with the inverse of an effective lensing efficiency /Jeff = piz) ■ 
P{z). Note that the source distribution p{z) has now to be 
renormalized over the range Zmin ^ z ^ 2max. 

In case of Euelid and all other ongoing and upcoming 
lensing surveys, photometric redshifts will also be available. 


Hence, we will only consider galaxies as sources for the lens¬ 
ing signal with redshifts z^in = 2phot > ^duster -1-0.15 where 
we choose an offset of 0.15 because the lensing contribution 
of sources close to the cluster redshift is negligible and for 
low redshifts the offset of 0.15 corresponds to the expected 
3a uncertainty in photometric redshift a^ = 0.05(1 -|- z). 

The tangential shear due to an NFW-profile can be ex¬ 
pressed as 

NFW/ N ENFw(a;) — E1 nFw(®) /ns 

7t (x) = -^ , (13) 

2-/crit 

i.e., as a scaled density contrast between the average surface 
density inside projected radius x and the surface density 
at radius x. However, observationally, it is only possible to 
measure the reduced tangential shear gr, i.e., 

9t = • (14) 

1 — K 

A parametrized model based on the equations above can 
then be used to derive the mass of the halo from the mea¬ 
sured shear signal. Here, we will also include the effect of cos¬ 
mic noise. This is important in order to derive more realistic 
uncertainties on the mass estimates (Hoekstra 2001, 2003; 
Dodelson 2004; Hoekstra et al. 2011). The most straight¬ 
forward implementation in order to achieve that is fitting 
the parametric model directly to a pixelized map of the two 
Cartesian projections of the reduced shear gx in the lens 
plane. These are related to gr through 

gi = gr cos{24>) (15) 

sin(2())). (16) 


For the implementation of cosmic noise contributions we fol¬ 
low Oguri et al. (2010) and calculate the 


2 

^ ^ [9a{0k) - ga{dk-,p)][C~%i3,ki 

q:,/3 = 1 k,l = l 

X [gfiW - > 


(17) 


where Greek indices run over the two components of the re¬ 
duced shear g and Roman indices run over the pixel positions 
{k,l = 1,..., Apixei). The matrix C denotes the covariance 
matrix and C~^ is its inverse. 

By minimizing the x^“value given the distortion ’data’, 
we find the best-fitting model parameters. Of course, with 
data we are referring to the shear components derived from 
a fiducial parametric model. 

We consider now two contributions in the covariance 
matrix: the dominating intrinsic ellipticity noise, i.e., shape 
noise, and cosmic noise due to large-scale structure along 
the line of sight. Thus, we write the covariance matrix as 

^ ^ ^shape ^LSS 


The intrinsic ellipticity noise between different galaxies is 
uncorrelated, thus, the shape noise covariance matrix con¬ 
sists only of diagonal terms 

= S^gSkl , (19) 

where 5ij denotes the Kronecker delta and 

Is the shape noise in pixel k. We estimate the number Nk of 

background galaxies per pixel as: 


Nk 



dzpBTc{z ), 


( 20 ) 
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where Apix stands for the area of one pixel and the source 
redshift distribution Psiciz) is the one given in equation 11 . 
We employ a number density of background sources of 
n = 30arcmin“^ which is expected for the Euclid survey 
(Laureijs et al. 2011). Furthermore, we assume an intrinsic 
ellipticity noise per galaxy of Uint = 0.25 per component. 

In contrast to that, large-scale structure along the line 
of sight introduces correlated noise which can be expressed 
as (Hoekstra 2003; Dodelson 2004) 

= \9k — 0i\) , ( 21 ) 

where denote the cosmic shear correlation functions. 
Since the universe is statistically isotropic, we express ^ as 
a function of the length r of the vector connecting the two 
positions Ok and 6i. In particular, the shear correlation func¬ 
tions constructed from the two shear components are given 
by 

^ii(r) = cos^( 2 <(>) 5 ++(r)-|-sin^( 20 ) 5 xx(r), ( 22 ) 

^ 22 {r) = Slid{2(1))^++{r) -|-cos^(2(())5xx(r), (23) 

62 (r) = ^ 2 i(r) = cos(2()))sin(2<^)[C+-H(r) -,fxx(j')] , (24) 

where (j) is the position angle between the coordinate x-axis 
and the vector r — Ok ^ Oi- ^++ and ^xx denote the tan¬ 
gential and cross-component shear correlation functions (e.g. 
Bartelmann & Schneider 2001), respectively. For the calcula¬ 
tion of these shear correlation functions we have to calculate 
the non-linear matter power spectrum Ps{k) folding in the 
source redshift distribution defined in equation 11. We cal¬ 
culate the non-linear matter power spectrum with the pub¬ 
licly available Boltzmann-code CLASS® (Bias, Lesgourgues & 
Tram 2011; Audren & Lesgourgues 2011). 

In all of our subsequent analyses, we £t the shear signal 
on a regular grid with constant side length A^pixei = 20 which 
is set to correspond to a square with side length 2 x 2 i?vir. 
Furthermore, we cut out the cluster centre in a square of side 
length 2 x 0.2 i?vir (we will refer to that more conveniently as 
fitting ’from 0.2J?vir to 2i?vir’ from now on). Larger scales 
than 2 are avoided since these are completely dominated 
by cosmic noise due to the smallness of the cluster signal 
there. In addition to that, two-halo contributions will also 
bias the NFW-fit on these scales (Becker & Kravtsov 2011). 
Smaller scales than 0.2f?vir are omitted since this is usu¬ 
ally done for practical purposes when dealing with real data 
In order to minimize the residual contamination by cluster 
members. Furthermore, the accuracy of shape measurements 
in high density and hence high shear regions is also an issue 
one tries to circumvent in practice. Finally, we also do ex¬ 
pect deviations from the simple NFW-profile on these scales 
due to effects of substructure (Becker & Kravtsov 2011). 

2.4 Results 

With the formalism outlined above we can now turn to de¬ 
termining the statistical precision of mass estimates from 
weak leasing measurements. Most importantly, we empha¬ 
size that realistic estimates of the statistical uncertainties 
must include cosmic noise (cf. equation 18). 

This is demonstrated by Fig. la and Fig. lb which show 

® Version 2.4.0 from www.class-code.net 


the likelihoods P{M) and P(c) as functions of the halo 
mass M and concentration c, respectively, resulting from 
the fitting procedure described above for three halo masses 
in the range 0.87 ^ M/(10 ^’’/i“^Mq) ^ 8.76, all at the 
same redshift 2 = 0.1875. In the derivation of P{M) we use 
the concentration-mass relation given in equation 8 . In the 
derivation of P(c) the mass was fixed at the fiducial value 
for each cluster and the fiducial concentration again derived 
through employing equation 8 . We then fit for the concen¬ 
tration c and demonstrate the importance of accounting for 
large-scale structure contributions in the data covariance. 

The dashed lines show the distributions only taking into 
account shape noise, whereas for the solid lines the effect 
of large-scale structure is additionally taken into account 
(cf. equation 18, 19, 21). This results in a significant broad¬ 
ening of the corresponding distributions P{M) and P(c), 
which in turn increases the uncertainties in the halo mass 
and concentration significantly. This effect is expected be¬ 
cause large-scale structure contributions reduce the weight 
of large scales on the weak leasing signal. Therefore, It is 
more pronounced for lower mass haloes since there the leas¬ 
ing signal is much weaker on larger scales compared to high 
mass haloes. 

We continue the analysis by looking further at the 
stacked weak lensing signal of clusters in eight mass bins 
from 6.78 x lO^®h“^M 0 to 2.70 x lO^®h“^M 0 and in four 
different redshift bins in the range 0 < « < 1.5. 

The relative error of the stacked lensing signal of a halo 
in each redshift bin i € 0, ...,3 and mass bin j € 0, ...,7 is 
estimated by 

^^teft/right ^ - Mj(2i)l , (25) 

l(AM)f+ AM5^“) 

"-“2 MAzi) ’ 

where the is calculated again as dehned In equation 17. 
The minimal is centred on the central value of each mass 
and redshift bin, respectively. 

In Fig. 2 we show the expected statistical relative un¬ 
certainty on the mass for a stacked weak lensing signal from 
galaxy clusters in each of the mass and redshift bins spec¬ 
ified above. The solid lines correspond to the total relative 
error including also large-scale structure contributions in the 
error budget. These are most important for clusters at low 
redshifts « ^ 0.75 and increase with increasing halo mass 
driven mainly by the lower abundances of these high mass 
clusters. 

In Table 1 we provide the quantitative values for the 
relative uncertainties in all mass bins for the lowest redshift 
bin since these ultimately set the requirements on the pre¬ 
cision. 


3 SYSTEMATIC ERRORS 

So far, our analysis did not include any sources of bias. How¬ 
ever, given the expected small statistical uncertainties as es¬ 
timated in the previous section (cf. Fig. 2 and Table 1), any 
source of systematlcs might become important. 

In general, the simplest parametrization of systematic 
deviations of the observed shear 7 obs from the true shear 7 
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M{h-^yiQ) 



(a) (b) 

Figure 1. (a) Likelihood as a function of halo mass M and concentrations derived through the concentration-mass relation of equation 8 
for single haloes of masses Mi = 0.87 x 10^^ h~^MQ (red), M 2 = 3.50 x 10^^ h~^MQ (blue) and M 3 = 0.87 x 10^^ (black). 

(b) Likelihood as a function of concentration c at fixed fiducial mass M for the same three haloes as in panel a) (fiducial concentrations 
are set by employing the concentration-mass relation of equation 8 . 

In both panels vertical lines indicate the fiducial values of mass and concentration, respectively. Dashed lines do not include, and solid lines 
do include large-scale structure contributions in the covariance (cf. equation 18). 



10 ^"* 10 ^® 


Figure 2. Relative uncertainties for masses of galaxy clusters as estimated from a stacked weak 
lensing signal for different mass and redshift bins. The solid lines take large-scale structure con¬ 
tributions into account whereas the dashed lines only include shot noise contributions (cf. Sec¬ 
tion 2.3). The eight mass bins (vertical, dashed lines) span a range from 6.78 x to 

2.70 X and the four redshift bins are in the range 0 < 2 : < 1.5. The horizontal grey 

line indicates the required upper bound on the multiplicative bias m for Euclid (cf. Section 3). 
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Table 1. Relative uncertainties in mass. 


Mio„(x1O14/i-1M0) 

MMgh(xlOl4/i-lM0) 

<^stat 

O'stat/'N/^ 

a 

CF/y/N 

N 

0.678 

1.075 

0.432 

0.0024 

0.605 

0.0034 

32509 

1.075 

1.703 

0.322 

0.0025 

0.477 

0.0037 

16861 

1.703 

2.699 

0.239 

0.0027 

0.378 

0.0042 

7999 

2.699 

4.278 

0.178 

0.0031 

0.302 

0.0052 

3363 

4.278 

6.780 

0.132 

0.0038 

0.242 

0.0070 

1200 

6.780 

10.75 

0.098 

0.0053 

0.196 

0.0106 

342 

1.075 

17.03 

0.073 

0.0086 

0.159 

0.0187 

72 

17.03 

26.99 

0.054 

0.0170 

0.130 

0.0410 

10 


The relative errors a/\/~N for each mass bin of Fig. 2 are given for the lowest redshift bin (i.e., 
•Smean = 0.1875). N refers to the number of haloes in each mass bin and Miow ^high denote 
the lower and upper value of each mass bin, respectively. 


can be written as (Heymans et al. 2006) 

7obs = (1 + rn )7 + c, (27) 

where m refers to the multiplicative bias and c to the ad¬ 
ditive bias (e.g. Huterer et al. 2006; Massey et al. 2013 and 
references therein). 

The multiplicative bias can arise from a variety of dif¬ 
ferent sources. For example, the shear is derived from an 
observed image which is a convolution of a point spread func¬ 
tion (PSF) of hnite size with the true shape of the observed 
object. This convolution introduces a multiplicative error. 
The additive bias, for example, can result from anisotropies 
of the PSF. We expect these to largely average out when 
referring to a stack of galaxy clusters which implies taking 
an average over many cluster-source pairs. 

Since the requirement on the multiplicative bias for Eu¬ 
clid is m < 2 X 10“® (Laureijs et al. 2011), we Hnd that the 
statistical uncertainties are still larger and thus this source 
of systematic can be neglected henceforth. Therefore, the 
primary observational sources of bias in weak lensing can 
both be neglected for our subsequent analysis in which we 
will focus then on other possible sources of bias instead. 

3.1 Photometric redshift bias 

For all ongoing and upcoming weak lensing surveys photo¬ 
metric redshifts of large samples of galaxies will be avail¬ 
able. Equations 10 and 12 show the dependence of the 
cluster lensing signal on redshifts of the cluster itself and 
all background sources or a source redshift distribution, re¬ 
spectively. Since it would simply require too much time to 
obtain spectroscopic redshifts for all galaxies in and around 
each cluster in these large samples of galaxy clusters, pho¬ 
tometric redshift estimates present the only feasible alter¬ 
native. However, the techniques for estimating photometric 
redshifts from multiband observations are not as precise as 
spectroscopic redshifts. They might even yield catastrophic 
outliers, for example, due to misinterpretations or confusion 
of emission or absorption features (e.g. Lya break; cf. Jouvel 
et al. 2011). 

This is also important for the weak lensing signal since 
cluster member galaxies that are scattered to higher red¬ 
shifts will now be treated as source background galaxies in 
the measurement of the weak lensing signal. In compari¬ 
son to the case of a perfect assignment of cluster members 


Table 2. The SN = 5 limiting (apparent) magnitudes for ex¬ 
tended objects used in the photometric redshift estimation. 


g r i 

z Y 

J 

H 

24.65 24.15 24.35 

23.95 23.3 

23.3 

23.3 


These limits assume DES ground based observations (g, r, i, z) 
and Euclid NIR observations (Y, J, H). The limiting magnitudes 
for extended objects are assumed to be 0.7 magnitudes shallower 
than for point sources. 

and background source galaxies the weak lensing signal will 
be diluted because erroneously scattered cluster members 
are not gravitationally sheared by the cluster. In this con¬ 
text another bias might arise due to cluster members be¬ 
ing intrinsically aligned although Sifon et al. (2015) showed 
no evidence for this. Note that any member galaxy scat¬ 
tered to a redshift in the foreground of the cluster does not 
change the lensing signal, though, because the photometric 
redshift distribution will be cut for all redshifts lower than 
^cut — .2^clus T 0.15. 

In order to quantify the level of contamination, we need 
an estimate of the number of cluster members we expect. For 
that reason, we adopt a Schechter-type cluster luminosity 
function (Schechter 1976) which follows the scaling relations 
with respect to cluster mass (or richness) derived from a 
large sample of groups and clusters from the Sloan Digital 
Sky Survey (SDSS) as presented In Hansen et al. (2009). For 
simplicity we only assume a passive redshift evolution of the 
luminosity function. The absolute magnitudes sampled from 
this cluster luminosity function are converted to apparent 
magnitudes which will then serve as input for an estimation 
of photometric redshifts for the cluster member galaxies. 

In general, observed magnitudes include noise from a 
variety of sources, for example, errors in the background 
subtraction, CCD readout noise, zodiacal light, and photo¬ 
metric errors. However, in our simulations it is sufficient to 
only use fixed magnitude limits. In Table 2 we show the 
signal-to-noise (SN) limits for different bands, combining 
DES optical photometry (e.g. Banerji et al. 2015) and Eu¬ 
clid near Infrared (NIR) observations (Laureijs et al. 2011). 
The photometric redshifts are estimated from the apparent 
magnitudes by employing a template based photo-z code de¬ 
rived from BPZ (Benitez 2000). A single redshift is estimated 
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by finding the peak of 




{h-Mz,T)y 


(28) 


where the sum is taken over the different filters i and differ¬ 
ent types of spectral energy distributions (SED) T for mim¬ 
icking different types of galaxies. Here /; is the observed flux 
in filter i, af. is the flux error, and fi{z,T) is a model flux 
constructed from the template library. In the estimation of 
the photometric redshifts, we do not include the Euclid vi¬ 
sual (VIS) band. Furthermore, no priors are included and 
the templates equal the ones used to generate the simula¬ 
tions. 

This yields a catalogue of photometric redshift esti¬ 
mates for cluster members for which we do know the true 
cluster redshift by construction. We cut the catalogue by 
imposing a detection limit for extended objects of < 

24.5 (lOa) in apparent magnitude following the requirement 
for lensing sources as given in Laureijs et al. (2011). After 
having applied this cut we use the catalogue to construct a 
photometric redshift distribution for the cluster members, 
4'cius(2), which we use to define the modified total source 
redshift distribution p(nod(^) Psr radial bin i; 

Pmod(«) OC ■ Psrciz) + iV*lus ' Pclua(z) . (29) 

Here denotes the total number of source galaxies per 
radial bin, i.e., = (1 —/!()?*)• Agfn. The number density 

of sources is again assumed to be n = 30arcmin“^ and with 
Afit we denote the area of an annulus defined by the borders 
of the radial bin. Note that the actual number of source 
galaxies per radial bin entering in the weak lensing analysis 
is again calculated as in equation 20 with psic replaced by a 
properly normalized Pmod • 

Moreover, we account for the fact that background 
source galaxies are also scattered into the cluster foreground 
by introducing /g))?*, the fraction of background sources scat¬ 
tered into the foreground of the cluster. These additional 
source scatterers intensify the effect of cluster galaxies scat¬ 
tered into the background twofold: firstly, the number of 
sources per bin is lowered which increases the uncertainty of 
the weak lensing signal and secondly the last term ofpl^adiz) 
modified by the scattered cluster galaxies gets a higher, rel¬ 
ative weight. We estimate from an additionally simu¬ 
lated source redshift distribution employing the photometric 
redshift estimation algorithm described above but using now 
luminosity functions from Marti et al. (2014). 

The source redshift distribution has the same functional 
form as in equation 11. We estimate the number of cluster 
members in the same radial bin i as 


Kius = fk-NX,{^2R,i,), (30) 

where is the fraction of cluster members in the radial bin 
i, i.e.. 


j«h.gh 


(31) 


The surface mass distribution 'E{R) is derived from the 
NFW-profile given in equation 5 (e.g. van der Burg et al. 
2014). The number of cluster members is estimated by in¬ 
tegrating the luminosity function up to the absolute magni- 



1.0 


z 


Figure 3. 

Upper panel: Effective lensing efficiencies /3efr(z) for the 
reference source background distribution (black, dashed line; 
cf. equation 11 ) and for a modified source redshift distribu¬ 
tion (red, solid line; cf. equation 29) for a cluster of mass 
M = 8.76 X 10^^ h~^ M 0 at redshift zi = 0.56 in the fitting range 
of the first radial bin, i.e., 0.2ilvir ^ r Si lilvir- The modified 
source redshift distribution consists only of elliptical galaxies 
(i.e., SED type “EllOl”). 

Lower panel: Relative difference between both effective lensing 
efficiencies. 


tude corresponding to the detection limit in visual apparent 
magnitude multiplied with the volume of a cylinder 

with radius 2Rvii and height 4i?vir. Note that the calcula¬ 
tions as presented above are slightly inconsistent with our 
approach of effectively fitting within squares enclosing the 
corresponding annulus. However, as can be seen in Fig. 4, 
the bias decreases with increasing radius and thus this ap¬ 
proach does not change our conclusions below. 

When we fit masses, we cut and normalize Pmod(z) over 
the range Zcut 5) z ^ Zmax, where Zmax should be formally 
set to infinity. We use here, however, a high enough redshift 
beyond which Pmod(z) = 0 Vz ^ Zmax- The fiducial model 
to which we compare in the fitting makes use of the same 
normalizations in Pmod(z), but with the cluster redshift dis¬ 
tribution Pciua(z) set to 0. 

We compare the effective lensing efficiencies (cf. equa¬ 
tion 12) derived from these two distributions over the fitting 
range of the first radial bin between 0.2 R^h ^ r ^ 1 i?vir 
in Fig. 3 for a cluster of mass M = 8.76 x 10^"^ h~^ Mq 
at redshift z = 0.56. The cluster redshift distribution is as¬ 
sumed to consist entirely of elliptical galaxies (i.e., SED type 
“EllOl”) in this case. In the lower panel we show the excess 
in effective lensing efficiency due to an excess in the clus¬ 
ter redshift distribution based on misidentified photomet¬ 
ric redshift estimates. In general, applying the foreground 
cut p(z ^ Zciuster + 0.15) = 0 already removes a large frac¬ 
tion of cluster members from the source background dis¬ 
tribution. However, there still remains a small fraction of 
cluster members with wrongly assigned (higher) photomet¬ 
ric redshifts which are then erroneously treated as members 
of the source redshift distribution. This affects the lensing 
efficiency as shown in the lower panel of Fig. 3. The mod- 
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ified lensing efficiency then leads to a bias in a weak leas¬ 
ing based mass derivation. Hence, we study the impact of 
this effect by repeating our previous weak lensing analysis 
by replacing the analytic expression from equation 11 for 
PsTc{z) with Pmod(2) for two fiducial clusters with masses of 
M = 0.88 X 10^'‘h-^Mg and M = 8.76 x 10^'‘Mg at 
11 equidistant redshifts between 0.1875 ^ «ciuster ^ 0.9375. 
We always include large-scale structure contributions in the 
error budget. Moreover, we consider two radial bins i per 
cluster, where the first bin spans radii in the range 0.2 i?vir ^ 
r ^ 1 iivir and the second bin extends from 1 f?vir to 2 Rvii- 
Typical numbers we find for the fraction of cluster mem¬ 
bers within these bins are ~ 50% and ~ 32%, respectively. 
The remainder of cluster galaxies is concentrated in the very 
centre of the cluster, i.e., between ~ 0 ^ r ^ 0.2 77vir. 

The results of this analysis are summarized in Fig. 4 and 
reveal that the impact of cluster members scattered into the 
lensing source sample is a minor concern for an individual 
cluster as indicated in the upper panels where we show the 
relative mass bias, b, as a function of true cluster redshift, 
2cius. When we compare the relative uncertainties for stacks 
of galaxy clusters, afy/N, as obtained in the previous section 
(cf. Fig. 2 and Table 1) with the level of expected relative 
bias b due to imperfect photometric redshift estimates of 
cluster members as shown in the lower panels of Fig. 4, we 
find that the bias is strongly dependent on the radial bin 
under consideration, as expected. For the first bin between 
0.2 i7vir and 1 i?vir the bias can be severe for certain combi¬ 
nations of SED type and redshift. In the second bin between 
1 Rvir to 2 Rvh the bias is always within the expected sta¬ 
tistical uncertainties though for all SED types and redshifts 
but it can amount up to 40% of the statistical uncertainties. 

The dependence of the bias on the SED type used for 
the galaxy templates in the derivation of the photometric 
redshift estimates is exaggerated for most SED types be¬ 
cause we always consider all members to consist of only 
one SED type. Galaxies with irregular SED types show the 
strongest bias whereas the lowest bias is found for elliptical 
SED types although it is then of equal strength as the sta¬ 
tistical uncertainties. In real clusters though, the majority 
of member galaxies will consist of (red) elliptical galaxies 
with rising fractions of (blue) spiral galaxies as a function of 
increasing cluster redshift (Butcher-Oemler effect. Butcher 
& Oemler 1984). Hence, the total bias will be dominated by 
a combination of the individual biases of these SED types, 
whereas contributions from irregular galaxies, which create 
the strongest bias in our analysis, will be much smaller in 
reality. Furthermore, a mass estimate derived over the full 
fitting range will be less affected by this bias since such an 
analysis corresponds to taking the average over both radial 
bins. The same conclusions will also hold for a measurement 
of cluster density profiles so that the photometric redshift 
bias will be negligible in this case, too. 

For high redshift clusters most of the photometric 
misidentifications will be scattered to redshifts in the fore¬ 
ground of the cluster, especially when considering the im¬ 
posed detection limits. Thus the ratio of the bias over the 
statistical uncertainties will flatten towards higher cluster 
redshifts (cf. lower panels of Fig. 4). 


3.2 Miscentring bias 


The position of the minimum of the cluster potential, its 
centre, is unknown in general. The position of the brightest 
cluster galaxy (BCG) or the centre of the X-ray emission 
can be used as a tracer for the cluster centre, but they do 
not have to coincide with the centre as determined by weak 
lensing, especially not in unrelaxed haloes. Hence, we need 
to consider a distribution in offsets. 

For a two-dimensional offset Vg we can calculate the 
azimuthally averaged convergence profile K.{r) (Yang et al. 
2006): 


1 r 


knfw 




r'^ + r‘^ + 2r ■ rs ■ cos 


(e)) 


d6», 


(32) 

where «nfw is the convergence for a spherically symmetric 
NFW-profile as given in equation 9. 

Based on SDSS-like mock catalogues for the maxBCG 
BCG finder algorithm, Johnston et al. (2007) found that 
the distribution of offsets follows a two-dimensional Gaus¬ 
sian distribution. Although the miscentring in these mock 
catalogues can only be caused by misidentifications due to 
the algorithm by construction, we still consider their results 
to be sufficiently good approximations since in practice one 
would also employ algorithms to determine miscentring off¬ 
sets in samples as large as expected to be found in the Eu¬ 
clid survey. Therefore, we follow Johnston et al. (2007) by 
assuming a two-dimensional Gaussian distribution of offsets, 

P{rs) = ^ exp(-0.5(rs/crs)^), (33) 


where the effective scale length was typically found to be 
Us = 0.42/i“^Mpc independent of cluster richness. The re¬ 
sulting convergence for miscentred clusters is then a convo¬ 
lution of the above equations which yields 

K‘{r)= J P(rs) «:(r|rs)drs . (34) 

From this convolved convergence we can derive the reduced 
shear for a miscentred cluster and employ the lensing 
formalism of Section 2.3. The total shear signal for a stack 
of miscentred clusters must be weighted by the fraction / of 
correctly centred clusters though, i.e., 

<?*°‘ = /-<; + (i-/)-/. (35) 


In Fig. 5 we show both the reduced shear signal and the 
convergence of a maximally miscentred halo (/ = 0 and 
as = 0.42h“^Mpc) and compare it to the expected signal 
of a perfectly centred halo of the same mass. The effect of 
miscentring is to dilute the shear signal in the inner regions 
of the cluster. Hence ignoring this will bias masses low with 
a dependence on the innermost radius used while fitting. 

Before considering strategies for mitigating the bias, we 
first examine its size further in Fig. 6a and Fig. 6b. Here, 
the miscentred shear signal is parametrized by the mass of 
the halo M and the two miscentring parameters. In Fig. 6a 
we fit a fiducial, centred cluster signal to the signal of a 
miscentred halo with a varying fraction of centred haloes 
(0 ^ 1) for three fiducial masses in the range 0.87 

M/(10^^ h“^Mg) ^ 8.76 (dashed lines). The width of the 
miscentring distribution is kept fixed at as = 0.42h“^Mpc. 
Miscentring introduces a bias in the recovered mass and in¬ 
creases with decreasing fraction of centred haloes /. Fur- 
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(a) Mtrue = 0.88 X IQi"* h-^ Mq 

0.2 Rvir r ss 1 Rvir 


(b) Mtrue = 0.88 X lO^^ h-1 Mq 

1 Rvir r s£ 2 iJvir 




(c) Mtrue = 8.76 X Mq 

0.2i7vir =£»■=£ lillvir 


(d) Mtrue = 8.76 X lO^^ h-1 Mq 

1 Rvir r- 2 Rvir 


Figure 4. 

Upper panels: Relative bias b due to imperfect assignment of photometric redshifts of cluster members which thus contaminate the source 
redshift distribution as a function of true cluster redshift ^duster^ SED type (i.e., ”E1101“ corresponding to elliptical, ”Sbc01“ and ”Scd01“ 
corresponding to spiral, and various ”Irr“ corresponding to irregular galaxy SED templates), and cluster mass. Panels (a) and (c) show 
the lowest and highest mass in the first radial bin, whereas panels (b) and (d) show the corresponding masses in the second radial bin. 
Lower panels: Comparison of the relative bias b from the upper panels with the corresponding relative uncertainties derived for a stack 
of galaxy clusters <j/y/N (cf. Fig. 2 and Table 1) again as a function of cluster redshift, SED type and cluster mass. Note the different scale 
of the y-axis between the radial bins. 


Table 3. Requirements on the precision of the fraction of centred haloes / and the relative error 
Acts/cts of the width as of the miscentring distribution. 


Mtrue (X 1014 /i-IMq) 

A/xnax 

(o-s = 0.42h-iMpc) 

V ■s / max 

{a-s = lh iMpc) 

\ ^ y max 

0.876 

0.0060 

0.019 

0.032 

3.500 

0.0171 

0.032 

0.036 

8.763 

0.0747 

0.088 

0.064 


We refer the reader to the text for details on the calculation of these requirements. 
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/ as{h ^Mpc) 


(a) (b) 

Figure 6. 

(a) Bias in halo mass introduced by miscentring as a function of the fraction of centred haloes / (for a fixed width as = 0.42/r“'^Mpc of 

the miscentring distribution) again with (solid lines) and without (dashed lines) large-scale structure contribution. The bias is shown for 
halo masses of Mi = 8.76 X 10^^ (blue), M 2 = 3.50 X 10^^ (red) and M 3 = 8.76 X 10^^ (black). 

(b) Bias in halo mass as a function of the width as of the miscentring distribution for a fixed fraction of centred haloes of / = 0.75 again 
with (solid lines) and without (dashed lines) large-scale structure contribution for the same halo masses as in (a). 



Figure 5. Reduced shear signals (blue) for a perfectly centred 
(/ = 1) halo (solid line) of mass M = 8.76 X 10^"* /r“^M 0 and its 
maximally miscentred (/ = 0 ) equivalent (dashed line) with as = 
0.42h“^Mpc (see text for details). The convergence is shown in 
red. The vertical dashed lines indicate the region which we include 
for fitting in our subsequent analysis (0.2i?vir ^ ^ ^ 2i7vir)- 


thermore, the bias is dependent on the fiducial halo mass 
and increases from high halo masses (of order 10^® 
to low halo masses (of order /i“^M©). The dependence 
of the bias on the fiducial halo mass (for fixed as) is caused 
by that the dilution of the shear signal due to miscentring 
is higher for low mass haloes. Taking then again large-scale 
structure into account increases the mass bias even further 
(blue, solid lines) because cosmic noise reduces the relative 
contribution of large scales to the shear signal. 

We repeat these calculations, but this time we keep the 
fraction of centred haloes fixed at / = 0.75 and leave the 


width of the miscentring distribution, as, free to vary. The 
results are shown in Fig. 6b for the same fiducial masses. The 
functional behaviour is this time more complex, but again it 
is apparent that including large-scale structure contributions 
(solid lines) also increases the bias caused by miscentring. 

Thus, the combined effect of large-scale structure and 
miscentring is to increase the uncertainties in the determina¬ 
tion of the halo mass and to introduce a mass bias (towards 
lower masses). By comparing the findings of Fig. 6a and Fig. 
6b to the relative uncertainties on the cluster mass as pro¬ 
vided in Table 1, we conclude that for the mass estimate 
of a single cluster the bias due to miscentring is negligible 
compared to the relative uncertainty a. For a stacked signal 
with statistics such as will be provided by the Euclid survey, 
this is no longer the case, since the statistical error a will 
be reduced by a factor 1/\/N (cf. Table 1). Hence, it will 
be necessary to determine the miscentring parameters / and 
as to within a few per cent in order to derive accurate mass 
determinations. 

In order to quantify this, we calculate upper limits for 
the precision to within which / and the relative error on 
as must be known such that the bias due to miscentring is 
smaller than the Icr-error in mass (cf. Table 1). The upper 
bound on the precision of / is thus defined as 

= . (36) 


Furthermore, the upper bound on the relative error Aas/as 
of the width as is given as 


Acts 

as 


max 


a f dAM \ 
'/N \ das 


(37) 


The values corresponding to the three fiducial masses Mtme 
of Fig. 6b including large-scale structure contributions are 
provided in Table 3. 
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Next, we explore if these conclusions also hold if we 
marginalize over the miscentring parameters expressing our 
lack of knowledge about them. Note as well, that we still 
assume equation 33 to hold which in itself is an important 
assumption in this regard. Hence, we explore the parameter 
space by fitting a miscentred shear signal to the shear signal 
of a fiducial model with mass M = 8.75 x / = 1 

at redshift a = 0.1875 for different values of M, Us, and /. 
All three parameters are varied for 50 different masses in 
the range 2.03 x 10^'‘ ^ M ^ 2.03 x 10^® /i"^M 0 , 40 

different values of (Js in the range 0 ^ CTs ^ 2 h“^Mpc, and 
20 different values of the fraction of centred haloes in the 
range 0 ^ ^ 1 , respectively. 

For the fitting we adopt again the formalism of Sec¬ 
tion 2.3 always including large-scale structure contributions. 
We have repeated these calculations for a second halo with 
fiducial mass M = 3.50 x 10^^ h“^M 0 (where we adjust the 
mass ranges accordingly). 

In Fig. 7a we compare the marginalized probability dis¬ 
tributions assuming flat priors on Os and / over the ranges 
indicated above for the two haloes with their perfect fit coun¬ 
terparts (cf. Fig. la). The marginalization over / and cjs, 
respectively, introduces a small bias due to the truncation 
of the prior on CTs at 2/i“^Mpc. However, the assumption of 
flat priors on both miscentring parameters is too pessimistic. 

Already with current data it is possible, particularly 
for the most massive clusters, to derive more realistic pri¬ 
ors for the miscentring parameters than fiat ones. Based on 
measurements of the displacement between the BCG and 
the maximum of the X-ray radiation from the intracluster 
medium from a sample of 53 massive galaxy clusters pre¬ 
sented in Bildfell et al. (2008), we have estimated errors on 
/ and (Ts by bootstrapping the distribution of measured off- 
centre radii. We find A/ = 0.04 and Acts = 0.01h“^Mpc 
independent of the number of bootstraps once this exceeds 
~ 1000. Interpreting these errors then as widths of Gaussian 
priors on / and cjs, respectively, yields qualitatively unbiased 
mass estimates when marginalizing over these priors instead 
of fiat priors, as shown in Fig. 7b. 

In the very near future a mapping of the X-ray sky will 
be carried out by eROSITA (Merloni et al. 2012) which is 
also aimed at detecting a large sample (~ 10 ®) of galaxy 
clusters. This sample can then be used to derive even better 
priors on the miscentring parameters which will mitigate the 
effect of the miscentring bias for Euclid weak lensing cluster 
masses entirely. Moreover, the combination of the eROSITA 
and Euclid cluster surveys will also allow to reduce the un¬ 
certainty on the hydrostatic mass bias between X-ray and 
weak lensing mass estimates so far that determining the sum 
of neutrino masses will be possible with cluster counts (cf. 
figure 12 from Planck Collaboration et al. 2015b). 


4 CONCLUSIONS 

We have investigated the level of statistical uncertainties and 
scrutinized various sources of systematic errors in the deter¬ 
mination of masses for stacks of galaxy clusters employing 
weak gravitational lensing. Throughout the analysis we have 
been focussing on future large area surveys, in particular the 
Euclid survey, so that the predominant source of statistical 


uncertainty, i.e., shape noise, is as small as possible when 
referring to (large) stacks of galaxy clusters. 

In Section 2 we have established the level of expected 
statistical uncertainties on the mass of (stacks of) galaxy 
clusters. We emphasize that contributions of cosmic noise 
must be included in realistic statistical uncertainties. How¬ 
ever, the level of expected statistical uncertainties is very 
low (cf. Table 1 and Fig. 2). Thus, at this level of expected 
statistical precision sources of systematic errors become im¬ 
portant which are still negligible compared to shape noise in 
current and ongoing surveys. Furthermore, the Euclid sur¬ 
vey poses an upper limit of m < 2 x 10“® (Laureijs et al. 
2011 ) on one of the predominant biases of weak lensing, i.e., 
the multiplicative bias m. The other predominant bias, i.e., 
the additive bias, is negligible when referring explicitly to 
stacks of clusters due to the averaging process involved in 
that kind of analysis. Hence, we expanded our analysis to 
other sources of bias instead. 

We have addressed one of these additional biases in Sec¬ 
tion 3.1 in which we investigated the effect on the accuracy 
of lensing masses when accounting for imperfect photometric 
redshift assignment to cluster member galaxies. Typically, a 
small fraction of these will be scattered to higher redshifts 
and thus mimic background source galaxies. This leads to a 
decreased weak lensing signal and thus requires the assign¬ 
ment of higher cluster masses during the fitting. We have 
shown that this bias is significant, especially for analyses 
using radii between 0.2 7?vir and IRvir- However, even in¬ 
cluding larger radii out to 2 Rvir will still require to prop¬ 
erly account for this bias in the full analysis and we strongly 
recommend to study the impact of this effect further in the 
context of more detailed simulations. 

Another source of systematic error is the effect of mis¬ 
centring, i.e., the ambiguity in the choice of an observational 
cluster centre. In general, a displacement of the observed 
cluster centre and the true cluster centre leads again to a di¬ 
lution of the shear signal and thus to higher mass estimates. 
Already with currently available data it is possible though 
to derive realistic priors for the miscentring parameters so 
that it will be possible to mitigate this bias entirely even 
In the case of Euclid when taking complementary missions 
such as eROSITA Into account. 

Finally, we want to emphasize that our analysis of these 
additional sources of systematic errors are all based on (sim¬ 
ple) analytic models. Eventually, these will have to be re¬ 
assessed by extensive numerical simulations in order to de¬ 
rive more realistic bounds and quantitative estimates. Our 
analysis, however, is meant to point to the relative impor¬ 
tance and order of magnitude predictions of these systematic 
errors in order to supply the community with a guideline for 
these future simulations. 
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(a) (b) 

Figure 7. Comparison of likelihoods for two haloes of masses M = 8.75 X 10^^ (left) and M = 3.50 X (right), 

respectively (blue, dashed lines) and likelihoods for the same mass haloes but marginalized over miscentring parameters ag and / (red, 
solid lines) assuming the following priors for both parameters: 

(a) Flat priors. 

(b) Gaussian priors. 

The grey, solid line shows the fiducial mass for each halo and we always account for large-scale structure contributions in the error budget. 
All haloes are at redshift z = 0.1875. 


matplotlib (Hunter 2007). Cosmology related calculations 
were performed using the Python package CosmoloPy^. 
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